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Abstract 

The propagation of electron beams carrying angular momentum in crystals is studied using a 
multislice approach for the model system Fe. It is found that the vortex beam is distorted strongly 
due to elastic scattering. Consequently, the expectation value of the angular momentum as well 
as the local vortex components change with the initial position of the vortex and the propagation 
depth, making numerical simulations indispensable when analyzing experiments. 
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I. INTRODUCTION 



Recently, electron beams carrying topological charge were created successfully in the 
transmission electron microscope (TEM)P21. The potential applications of these vortex 
beams (as they were called) seem endless — the speculation is that they will become as ver- 
satile a tool as vortices in opticd^. In particular, proposed applications range from strongly 
enhanced energy-loss magnetic chiral dichroism signal over the probing of chiral structures 
to the manipulation of nanoparticles^. 

One essential aspect of electron beams has been overlooked so far, though, which is of 
paramount importance to real experiments. It is well-known that — due to the strong and 
long-ranged Coulomb interaction — the probe electrons in a sample undergo strong elastic 
scattering which depends heavily on aspects such as beam orientation or sample thicknesd^. 
This can be exploited, as in high-resolution TEM (HRTEM)f^, atom location by channeling- 
enhanced microanalysis (ALCHEMI)^, or energy-loss by channeled electrons (ELCE)^^. 
Recently, the new RSTEM technique was used to map the effect of elastic scattering on a 
focused electron beam with sub- Angstrom precisioii^. 

In other cases, such as energy-loss magnetic chiral dichroism (EMCD)-^^, it is usually 
detrimental to the signal strengtlP^. In all cases, however, it plays an important role, and 
the same is to be expected for electron vortex beamd^. 

II. SIMULATIONS 

In this work, we concentrate on the propagation of electron vortices that have been 
focused on the specimen, i.e., vortices produced as a convergent beam hitting the specimen. 
This corresponds to the situation of a holographic or phase plate in the condenser system 
of a TEIVpSHSl^ ii must be emphasized, however, that the same results also apply to vortex 
beams that emanate from some point within the sample, e.g., after a chiral excitation as in 
EMCD. 

In the case of a vortex produced by a condenser aperture, we can model the wave in the 
condenser plane a^ 

^(g,0) = e^-^e(w-g), (1) 

where m = 1 is the topological charge, (f) is the polar angle in the condenser plane, q is 
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FIG. 1. Projection of an Fe crystal in [0 1] direction. Atoms in different layers are drawn as disks 
and squares, respectively. The positions mentioned in the text are marked by capital letters. The 
edge length of the square is the lattice constant of Fe, namely 2.859 A. 

the radial vector in the condenser plane, g^ax is the radius of the aperture (which defines 
the maximum convergence angle), and Q{x) is the Heaviside function. In the focused case 
(ignoring aberrations) , the wave function ip hitting the sample is given by the Fourier trans- 
formation of t/S, which corresponds to the Hankel transformation in polar coordinated^ 
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^(^,^) = Ji{qr)qdq 

i<?maxe ^j^(^q^^^r)Ho{qm^^r) - Jo{qm^^r)Hi{q^^^r)) 



(2) 
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Here, we used the Bessel functions of first kind (J^) and the Struve functions (i^m)- 

As model system, we have chosen Fe. On the one hand, it has a simple crystal structure 
(body centered cubic), on the other hand it is an easily accessible sample material and of 
general interest as a ferromagnetic model system?^. In particular, we investigated the [0 1] 
zone axis which is a high-symmetry zone axis as can be seen from fig. [T] 

For the simulations, the focused vortex beam was positioned on four different positions on 
the crystal (see Fig. [T]): on an atomic column (A), between the two closest atomic columns 
(B), as far away from atomic columns as possible (C), and in another, non-symmetric place 
(D). In aU cases, a thickness of 200 A and a spherical aberration coefficient Cs = were 
assumed. Furthermore, an energy of 200 keV, an angular momentum of h and a conver- 
gence semi-angle of 15 mrad (corresponding to a waist radius of 0.67 A) were used for the 
incident electron beam. 

This electron beam was subsequently propagated through the crystal using the multislice 
approachP^. Fig.[2]shows the propagated wave functions at the exit plane. From them, it is 
obvious that (a) the wave functions do no longer show the typical vortex phase and amplitude 
structure^^ in general, (b) the structure — and in particular the phase distribution — depends 
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FIG. 2. Wave functions at the exit plane of the 20 nm thick Fe sample for the positions A, B, C, 
and D as defined in the text and in fig. [ij Column 1 shows the real-space amplitudes (brightness) 
and phases (hue, as specified in the color wheel on the right), column 2 shows the corresponding 
intensities, column 3 shows the reciprocal-space amplitudes and phases, and column 4 shows the 
diffraction intensities. The 1 and 2 show a (10 A)^ subset of the simulated data, while the columns 
3 and 4 show a (2.45 /A)^ subset. 



critically on the initial position of the focused vortex with respect to the crystal, and (c) that 
we generally see curved structures as one would expect for vortices, but which are absent in 
normal STEM12o1^ 

In addition to the curved features between the atomic columns, the simulations shown 
in fig. [2] also reproduce the channeling and dechanneling effects usually seen with focused 
electron beams^^. As a consequence, one sees significant intensity not only close to the 
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initially illuminated atomic column, but also on the nearest and next-nearest neighboring 
columns. Moreover, the central vortex phase structure is relatively well-preserved in-between 
the atomic columns, even after a 200 A thick sample. In close proximity of the atomic 
columns it is heavily distorted, however, which also leads to very different local phases at 
and around neighboring atoms. 

It must be noted here that all the images in fig. [2] display the wave functions directly 
at the exit plane of the sample. If one wants to observe such patterns experimentally in a 
microscope, one has to use some post-specimen lens system to transfer those wave functions 
to the screen. These lens systems are usually not ideal and hence can introduce artifacts 
or obscure some details. Since the important object for all interactions in the sample is the 
wave function, however, it is more appropriate to show this directly, rather than the image 
one would get on a camera. 

Recently, the feasibility of recording images similar to those shown in fig. [2] was also 
demonstrated experimentally using the new RSTEM techniqu^^^. In combination with a 
next-generation, double-corrected TEM, this technique can be used to investigate the real- 
space image of a focused electron beam with sub- Angstrom resolution after it has traversed 
the crystal. 

III. ANGULAR MOMENTUM 

One of the most important characteristics of a vortex beam is angular momentum. As 
helical wave functions are eigenstates of the operatoi'^, they have a well-defined, integer 
angular momentum quantum number m which is a constant of motion in free space. In fact, 
it is precisely this quantum number that defines the "vorticity" of the beam. But how does 
it behave inside the sample? 

Fig. [i] shows the evolution of the expectation value (L^) of the angular momentum as a 
function of depth z in the crystal. It is clearly obvious that (a) it changes significantly as 
the probe propagates through the crystal, (b) that it is generally not an integer and can 
even surpass the initial value of /i, and (c) that it is position-dependent. 

All of this can be explained by looking at the properties of L^. While this operator 
commutes with spherical operators such as the free-space Hamiltonian, it does not gener- 
ally commute with the crystal Hamiltonian in which the spherical symmetry is broken by 
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FIG. 3. Expectation value (L^) in Fe [0 1] for the four positions given in the text and in fig. [l]as 
a function of penetration depth z. The curves for the positions B and D overlap. 

the crystal potential. Thus, (L^) is constant in free space, but can and will change inside a 
crystal. This also implies that the probe inside the crystal will evolve into a (coherent) super- 
position of different eigenstates of over time, accounting for the non-integer expectation 
values. 

Physically, this can be understood as an exchange of angular momentum between the 
probe beam and the crystal. Because of the very different masses of probe electron and 
sample as well as the oscillatory natur^^of (L^) in fig-jsj the change of the angular momen- 
tum of the sample usually is not measurable. For nanoparticles of suitable sizes, however, 
this could produce appreciable angular velocities. 

From this physical interpretation, the position dependence is understandable as well. If 
the (Angstrom sized) probe is positioned directly on an atom column, it usually channels 
down along that column, whereas a similar probe positioned between columns will first 
spread until it hits several columns at once and will subsequently channel down primarily 
along those^^. Therefore, the amount the focused probe beam interacts with the sample 
obviously depends on its position, as is readily visible in fig. [3| 

IV. VORTEX DECOMPOSITION 

In the previous section, we have seen that inside the sample, a pure vortex state evolves 
into a coherent superposition of such states. In fact, the occurrence of different vortex 
components is not an effect unique to a crystal — even in free space, a vortex wave can be 
described by a single exp(im(/;) function (with polar angle (/;) only in a single cylindrical 
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FIG. 4. Expansion of an m = 1 vortex incident on position A for the positions A-D given in the 
text. Only the most important components ml — —1, 0, 1, 2 are plotted. The scale is the same in 
all subfigures. The vertical line marks the nominal beam waist radius that was used to calculate 
the depth plots in fig. [5j 



coordinate system, namely the one for which the z axis coincides with the propagation axis 
of the phase singularity in the center of the vortex. If the vortex is displaced relative to the 
origin of the coordinate system, again several components come into play, as can readily be 
seen from fig. [4j Only if we describe the m = 1 vortex centered on position A around that 
very position do we retrieve exclusively an m = 1 component. In all other cases, additional 
components are mixed in, including m = Airy-like components and even opposite m = —1 
components. 

In free space, all these components add up to a single (displaced) vortex wave which has a 
well-defined and — most importantly — constant angular momentum!^. In a crystal, however, 
this is no longer the case, and consequently the relative weights of the different components 
change with the penetration depth in the crystal, just as (L^) changes. 

Fig. [5] shows how the most important local amplitude components around the atom at 
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FIG. 5. Dependence of the local amplitudes of the most important vortex components m = 
— 1,0,1,2 for vortices positioned at the positions A-D given in the text on the depth z in the 
crystal. The local amplitudes were calculated at a radius of r = 1.08 A (which corresponds to the 
nominal beam waste radius) around the atom at position A. The scale in all subfigures is the same 
as in fig. [i) 



position A change with depth, dependent on the vortex position. For this, the vortex was 
propagated with an initial position A-D. At each depth, the wave function was decomposed 
in vortex components around the atom A, and the amplitude at the radius r = 1.08 A 
(which corresponds to the nominal beam waist radius) was extracted. As such, z = in 
fig. [5] corresponds to r = 1.08 A in fig. jij 

It is readily visible from fig.[5]that the local vortex amplitudes change with thickness, but 
the precise behavior depends crucially on the initial vortex position and the vortex order 
of interest. The oscillatory behavior seen best in the curves for A/m = l,D/m = l, 
and D / m = 2 is reminiscent of the pendellosung of conventional Bloch waves, although an 
analytical equation for its description is unfeasible because each vortex contains infinitely 
many Bloch wave components. 
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V. CONCLUSION 



The calculations shown in the previous sections show that one has to be careful when 
predicting or interpreting the effects of vortices if they pass through matter. In particular, 
statements one usually takes for granted when dealing with vortices — such as the conserva- 
tion of topological charge or equivalent ly angular momentum — are only valid in free space. 
Inside a sample, electron vortex beams can be heavily distorted, which will change the com- 
ponents present in a local environment around the atom. As a consequence, one can have 
m = and m = 2 beams in the vicinity of a scattering atom, even if the original vortex had 
m = 1. 

This creates new challenges for the use of vortices in, e.g., EMCD, but at the same time 
provides new exciting possibilities. These range from the conversion of vortex beams from 
one topological charge to another to the transfer of angular momentum to nanoparticles. 
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